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Abstract 

An analog of the lattice KdV equation of Nijhoff et al. is constructed on a hexagonal lattice. The 
resulting system of difference equations exhibits soliton solutions with interesting local structure: 
there is a nontrivial phase shift on moving between adjacent lattice sites, with the magnitude of 
the shift tending to zero in the continuum limit. 

o 

^ Bobenko, Hoffmann and Suris BHS have recently proposed that some effort should be made to 
^understand integrable systems of difference equations on lattices other than 7? . In addition to being of 

^theoretical interest, it is possible that in certain cases this could have physical significance MER,^KMG. 
'^This letter examines the analog of Nijhoff et al.^s lattice KdV equation NC^me on a hexagonal lattice, 

j>which we call "hexaKdV" for short. The soliton solutions of hexaKdV are constructed. HexaKdV 
^solitons undergo a nontrivial phase shift on moving between adjacent lattice sites, with the magnitude 

C^o'f the shift tending to zero in the continuum limit. 

Nijhoff et a/.'s lattice KdV equation involves a single field b defined at the vertices of the standard 
lattice Z^. The equation gives a relation between the values of the field at the 4 vertices of each 
fundamental lattice plaquette. Writing bi for 6„,m, 62 for bn+i,m, &3 for bn+i,m+i and 64 for 6„,m+i (so 
bi,b2,bs, 64 are values of the field as we go around a fundamental plaquette), the equation takes the form 
rne 

\ &1&2 + O2O3 - O3O4 + O4O1 = . (1) 

Two features of this equation are important in the sequel. The first is that the equation is unchanged 
on replacing bi by 62, ^2 by 63, 63 by 64, 64 by bi, hhj k and k by —h. This symmetry is associated with 
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rotating the lattice through 90°. Indeed, the equation can be written in the evidently symmetric, and 
rather more compact, form 



where R is the rotation (or replacement) operator. The second important feature to note about standard 
lattice KdV is that the number of variables and the number of equations is properly balanced. Each 
equation relates between 4 different values of the field, but each value of the field appears in 4 different 
equations (as each vertex of the lattice belongs to 4 fundamental plaquettes). Thus lattice KdV is a set 
of equations that can be expected to yield a solution, given suitable initial/boundary data. 

The hexaKdV system which will shortly be constructed also displays a rotation symmetry. Consistent 
with this, each equation in hexaKdV relates the 6 values of the field at the vertices of each fundamental 
hexagonal plaquette. But now a problem arises: Each vertex of the lattice evidently belongs to only 3 
fundamental plaquettes, so it seems there are twice as many variables as equations (each equation calls 
6 variables, but each variable only appears in 3 equations). The simplest imaginable solution to this 
problem is that there should be two equations associated with each plaquette. In fact there is a system 
of four equations associated with each plaquette, but they are degenerate and should be considered like 
two. 

The hexaKdV system will be constructed on the hexagonal lattice with vertices 



where = e^*'^/^, and g, h, k are arbitrary positive reals. See figure 1. The vertex riiq + n2h0 + n^k0^ 
will be referred to in the sequel simply as "the vertex ^1,^2, 713" . There are 3 kinds of edges, parallel to 
1, 0, 0^ respectively. Following the ideas of BHS and NC, we look for a G'L(2)-valued function ^ni,n2,n3(0 
defined on the vertices of the lattice and dependent on the spectral parameter A, satisfying the following 
equations, one for each edge in the lattice: 



The consistency condition arising from the two different ways to go around a fundamental plaquette 
(rii, 712, n-i, ^2 + 1, 713 ^ rii, 7^2 + 1, 713 - 1 ^ rii + 1, 7^2 + 1, 713 - 1 or 7ii, 712, 713 ^ 7ii + 1, 712, 




(2) 



{niq + n2h0 + n3k0'^ : rzi, rz2, e Z, rzi + 712 + ris = or 1 } , 



(3) 
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Figure 1: The hexagonal lattice 
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Figure 2: Two ways around a fundamental plaquette 



rii + 1, n2, ^3 — 1 ^ ni + 1, 712 + 1, ris — 1, see figure 2) can be written 



1-qh q ' 

^ql + b-s- h- qhzhi 1 + 563 , 



kl + b2 - - A;62&3 I + kb2 , 



1-/162 h 

hi + bi - b2 - hbib2 1 + hbi 



(5) 



1 — hb^ h 
/i/ + 65 — 64 — hb5b4 1 + hb^ 



1 - qba 
^ql + bi - be - qbibe 



q 

1 + qbi, 



1 — kb^ k 
kl + be — b^ — kb^b^ 1 + kb^ 
where here for brevity bi stands for 6„i,„2,„3, 62 for 6„i,n2+i,n3) ^3 for &ni,n2+i,n3-i> ^4 for 6„,+i,„2+i,„3_i, 65 
for 6„j+i „2,n3-i and 65 for &ni+i,n2,n3) these being the 6 values of b at vertices of a fundamental plaquette. 

A symbolic manipulator was used to multiply out cons. Four equations emerge. Since full rotational 
symmetry has been broken by the choice of a starting and ending point in the fundamental plaquette, 
the equations are not rotation invariant, but this can be rectified by taking suitable linear combinations 
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and just the final symmetric form of the equations will be presented here. Writing R for the 60°rotation 
operator, implemented by replacing bi by 62, ^2 by 63, 63 by 64, 64 by 65, 65 by be, b^ by 61, h hy k and 
A; by g and q by —h, the equations are: 



(6) 



(l + R + R'' + R' + R^ + R') [i^l— M! _ 6^62(&i - &2)) = , (7) 



(/ _ + _ ^3 ^ ^4 _ ^5 j (^^^^^ - 6l64(&2 + h) (8) 

, 6163 + 6164 + b2b3 - 62^5 - hh - &3^'5 



h 



= 0, 



(l-R + R'^-R^ + R''- i?^) (|^^^^^^^ + ^1^2(^1 + &2 - ^1^2 - 263^4 - hh - 6265) (9) 

bl-bl + 6i(26i + 65) - 6^(64 + 262) + 261(64 - ^2) + 2(6162^6 + fe2&3&5 - fci&2&3 - bib^bM 

h 



. 



Note the similarity of the first equation, hkdvl, with the standard lattice KdV equation dkdv. An 
analysis of the hexaKdV system hkdvl-hkdv4 is quite straightforward with a symbolic manipulator. 
The first equation, hkdvl, is linear in 61, and can be used to eliminate 61 provided its coefficient, 
be — b2 + — ^, does not vanish. Proceeding along these lines it can be shown that the solution set 
consists of two components, characterized by whether or not &6 ~ ^2 + ^ — ^ vanishes. One component 
is a dimension 3 set of solutions given by solutions of the 3 linear equations 

66-&2 + i-i = 

fe2-fe4 + ^ + i = } . (10) 

(^-^)6i + (^-^)63+(^-^)5, = (i-i)(iH-i)(iH-i) 
The other component is a dimension 4 set of solutions in which 61 is determined by 



63 64 - 63 62 - 65 64 + b,b, + ^^^^ + -''3+^fe4-&6 + 



^1 = 1 l—l 1 ' (11) 



h,-b2 + l 



and 62, &3, &4, &5, be satisfy the single constraint 

= b2b^be + be 63 64 - 62 65 64 - be 63 62 + 65 64^ - 63 64^ + 63 64 62 - 65 64 be 

^ 64 ^6 - - ^6 62 + 65 be - 63 66 + 64 62 + 63 64 - 65 64 _^ 65 + 65 - ^3 - 2 64 + 62 

a hk 
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^ -63 65 + 65 &6 - 2 65 64 + 2 63 64 -hb2 + 62 ^5 , -^2 -h + be + b^ 

^ 63 ^4 + ^2 ^5 + ^4^ - &5 ^4 - ^4 ^2 - ^3 ^2 - ^4 ^6 + ^2 _^ -^e + 2 64 - 62 - ^3 + &5 

^^+hr,-W, + b2 



A;2 

It can be checked that the 4 dimensional set of solutions includes the 3 dimensional set obtained by 
rotating the solution of 3dset, i.e. the solutions of 



(13) 



Thus, to summarize, we have arrived at the hexaKdV system. This consists of the system hkdvl- 
hkdv4 of 4 polynomial equations in 6 variables on each plaquette of the lattice. Because the system 
hkdvl-hkdv4 is degenerate, in the sense that it has a 4 dimensional solution set, it is more appropriate 
to think of hexaKdV as specifying 2 constraints for each plaquette. Thus, as explained above, there is 
a proper balance between the number of variables and the number of constraints in hexaKdV. 

Soliton solutions. A direct computation shows that hkdvl-likdv4 has the following two parameter family 
of solutions: 

bi — Ctanh(2;) , 

62 = Ctanh(2 + tanh-^(/iC)) , 

63 = Ctanh(2 + tanh"^(/iC) +tanh-^(A;C)) , 

64 = Ctanh(z + tanh-^(/iC) +tanh-^(A;C) +tanh"^(gC)) , 

65 = Ctanh(^ + tanh-^(A;C) +tanh-^(gC)) , 
be = Ctanh (2; + tanh"^(gC)) . 



Here C, z are arbitrary. Such solutions on an individual plaquette can be pasted together to give a full 
solution of hexaKdV of form 

bni,n2,n3 = C tanh (^ni tanh^"*^ (gC) + n2 tanh^^(/iC) — n3 tanh^"'^(A;C) + , (15) 

where again C, z are arbitrary constants. To get some understanding into the nature of the soliton 
solution, it is necessary to write it as a function of the standard Cartesian coordinates x and y of the 
vertex ni, ^2, 713. These are given by 

1 ^/3 
X ^ mq - -{n2h + nsk) , y ^ -^{n2h - Usk) . (16) 
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Writing ni,n2, in terms oi x,y and the quantity 



s = ni+n2 + n3 (17) 

(which is or 1), the sohton becomes 

b{x, y, s) = Ctanh(i:>(x + cy) + Es + z) , (18) 

where 

^ _ {k + h)teinh-\qC) + htanh-\kC)-kteinh-^hC 

qk + qh + hk ' 
{h- k) ia.n\r\qC) + {h + 2q) tanh"^ (/eC) + {k + 2q) tanh"^ hC 
^ ~ y^(^{k + h)tanh-\qC) + htanh-\kC)-ktanh-^hC) 

^ _ kh taiuh'^ (qC) - qh tanh~^ (kC) + kq taiuh'^ hC 

~ qk + qh + hk ' ^ ^ 

Before discussing this result, note that for sufficiently small q, h, k (for which tanh~^(gC) ^ qC etc) 

D ^ ^(^±^, (22) 
qk + qh + hk ^ ' 

^ _ qk + 3qh + 2kh 

V^q{k + h) ' ^ ^ 

E ^ '-^ . (24) 

qk + qh + hk ^ ^ 

In the form sol2, the meaning of the sohton solution is quite clear, c is the speed of the soliton, C 
its amplitude, and equation sa expresses the speed-amplitude relation. The dependence of D on C, 
through wa, shows that the width of the soliton also depends on the speed (or amplitude), a common 
phenomenon seen, for example, in the KdV equation. The novel feature of hexaKdV solitons is, however, 

the dependence of sol2 on s. This means that there is a phase shift between the soliton solution on "even" 
(s = 0) and "odd" (s = 1) sites of the lattice. This interesting phenomenon distinguishes hexaKdV from 
standard lattice KdV. Equation hm3 shows that as the continuum hmit is approached (g, h, k tending 
to 0), the dependence on s becomes very weak, and indeed vanishes in the continuum. Equations liml 
and lim2 show that as g, h, k tend to with constant ratios, the speed and width of the solitons tend to 
constant values, independent of amplitude. 

The continuum limit. Is hexaKdV in any sense a discretization of a PDE? Consider making the following 
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replacements in hkdvl-hkdv4: 



h{x-\q + \h-\k,y-^h-^k) 
hlx-\q-\h-\k,y + ^h-^k 



X 



hix + \q 



\h 



4 



ik,y 



lk,y + ^h + fk 



4 



^h+^k 



(25) 



b[x + lq+ \h - \k,y - ^h - ^k) . 



Expanding the 4 resulting equations in a Taylor series in q, h, k and retaining only leading orders (as- 
suming /i, g, A; all to be of the same order) gives only 2 distinct equations. Equations hkdv3 and likdv4 
both give 



qk + 3qh + 2kh 



bx 



(26) 



^ VSq{k + h) 

This is consistent with the results on soliton solutions, as for small q,h,k the soliton speed is given 
by c as in lim2, and this is precisely the factor that has just appeared in weq. However, both hkdvl 
and hkdv2 reduce to another PDE, which is only consistent with weq ii q — h. So far we have no 
understanding of why it should be necessary to impose such a constraint for a consistent continuous 
limit, or whether there is maybe some reason to ignore it. 



Dual Discretizations of PDEs. Although it is a digression from the main topic, since it has been shown 
that hexaKdV is a discretization of weq, at least in the case h — q, we briefly raise the question of 
how ab initio one might go about discretizing this (or another) PDE on a hexagonal lattice, with the 
intention of getting equations relating the 6 values of the field around each hexagonal plaquette. 

The standard finite difference approach for discretizing a PDE involves writing a discrete equation to 
approximate the PDE at each vertex of the lattice being used. The discrete equation involves the values 
of the functions appearing in the PDE at the relevant vertex, as well as the values at neighboring vertices. 
The balance between the number of equations obtained and the number of variables is automatic, as 
both equal the number of lattice vertices in the relevant domain. 

Here we follow a different "dual" approach. The aim is to approximate the PDE on a fundamental 
lattice plaquette (or, more precisely, at the center of gravity of a lattice plaquette), using the values 
of the function on lattice vertices. In the case of a hexagonal lattice, this means it is necessary to 
approximate derivatives of b at the point {x, y) using the 6 values of b that appeared in fepls. This is 
straightforward. For example it can be shown that under the replacements fepls 

^ b^{x,y) + 0{q',h',e,hk,qk,qh), (27) 
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Thus the PDE by — Mb^/y/S can apparently be discretized by the difference equation 

65 - 62 - &6 + ^3 _^ 64 - &i - 65 + h ^ ^ b4-bi + be- 63 ^^g^ 



k h q 

However, there is now a counting problem of exactly the type mentioned before. The PDE has been 
replaced by a single difference equation for each plaquette, with the difference equation "calling" 6 values 
of the field. But, each value of the field only appears in 3 equations, so the number of variables is twice 
the number of equations. The resolution of this is that there is a constraint that needs to be imposed. 
Working to the same order as the approximations fdl-fd2 it can be checked that 



Thus equation papl on each plaquette must be supplemented by the constraint 



(30) 



1 + -] (64 - 61) -d + l) ib, - 62) + f - + 7 1 (^6 - &3) = . (31) 



h qj \k hj \q k^ 

This resolves the counting problem. 

Dual discretizations of PDEs may well be appropriate in a variety of settings. It is certainly no 
surprise that they arise in discretizations of equations in the KdV hierarchy, which are, in a natural 
way, zero curvature equations. In conclusion of this section, note that the constraint equation pap2 
looks very similar to what is obtained by applying all the necessary rotation operators just to the term 
(61 — b2)/h in hkdvl. But the differences turn out to be significant; for hexaKdV there does not seem 
to be a natural division of the equations into a constraint part and a dynamic part, as there is with the 
simple discretization comprised of papl and pap2. 



OctaKdV and other lattices. The work presented in this letter can be extended to other lattices. We 
report briefly on the extension to a lattice of alternating octahedra and rectangles (with the octahedra 
having 4 pairs of parallel sides of equal length). On the rectangular plaquettes standard lattice KdV 
must hold. A simple counting argument then shows that there should be 3 equations (in 8 variables) 
associated to each octagonal plaquette. Looking at the analog of cons with 4 matrices in each side gives 
a complicated system of 6 polynomial equations in 8 variables, with the first being an obvious extension 
of dkdv and hkdvl. We hypothesize that this system has a 5 parameter solution set, and thus should be 
thought of as just 3 constraints. As of yet we have not been able to verify this in general, the memory 
requirements for the algebra being substantial. We have, however, looked at a number of reductions, in 
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which one variable is assumed to vanish; in such reductions there are 4 parameter solution sets, which 
is consistent with the hypothesis. 

It may be possible to obtain some general results on the algebraic system obtained from the gen- 
eralization of cons with n matrices on each side, relevant to a plaquette with 2n sides. Presumably 
soliton solutions on more general lattices will also exhibit "small scale" phase shifts as we move round 
a plaquette, as there are for hexaKdV. These phase shifts may well be indicative of lattice structure. 

It remains a mystery whether there is a way to formulate integrable systems of difference equations 
on lattices with plaquettes with odd numbers of sides, particularly triangular lattices. 
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